Nonlinear resonant wave interaction in vacuum 



Gert Brodin, Daniel Eriksson 

Department of Plasma Physics, Umed University, 
SE-901 87 Umed, Sweden 

Mattias Marklund 
Department of Electromagnetics, Chalmers University of Technology, 
SE~^12 96 Goteborg, Sweden 

The basic equations governing propagation of electromagnetic and grav- 
itational waves in vacuum are nonlinear. As a consequence photon-photon 
interaction as well as photon-graviton interaction can take place without a 
medium. However, resonant interaction between less than four waves can- 
not occur in vacuum, unless the interaction takes place in a bounded region, 
such as a cavity or a waveguide. Recent results concerning resonant wave 
interaction in bounded vacuum regions are reviewed and extended. 
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I. INTRODUCTION 

Most examples of nonlinear wave phenomena occur as a result of some nonlinear property 
of a medium. In particular this applies to electromagnetism since Maxwell's equations are 
linear. Still, nonlinear interaction between photons in vacuum may occur as a result of scat- 
tering processes involving vacuum fluctuations, as described by quantum electrodynamics 
(QED). An effective theory for this can be formulated in terms of the Euler-Heisenberg La- 
grangian [1,2]. It should be noted, however, that resonant interaction between less than four 
waves requires that the waves are parallel, in which case the nonlinear QED coupling van- 
ishes [3,4]. Similarly, gravitons and photons couple in vacuum, as described by the Einstein- 
Maxwell system of equations. But the dispersion relation implies that the waves must be 
parallel, in case three wave interaction should be resonant. Just as for photon-photon inter- 
action in vacuum, however, this condition means that the wave-coupling vanishes [5]. On 
the other hand, the situation is changed if some or all of the interacting waves are confined in 
bounded regions such as waveguides or cavities [5-7] . Below we will review and extend recent 
results concerning resonant photon-photon scattering as well as photon-graviton scattering 
in bounded regions. 



II. PHOTON-PHOTON SCATTERING 

According to QED, the non-classical phenomenon of photon-photon scattering can take 
place due to the exchange of virtual electron-positron pairs. This give rise to vacuum 
polarization and magnetization currents, and an effective field theory can be formulated in 
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terms of the Euler-Heisenberg Lagrangian density [1,2] 



^EH = £o^ + £(4F 2 + 7£f 2 ) , (1) 

where £ = 2a 2 elh 3 /45m^,c 5 , & = \(E 2 — c 2 B 2 ), = cE • B, e is the electron charge, 
c the velocity of light, 2nh the Planck constant and m e the electron mass. The latter 
terms in (1) represent the effects of vacuum polarization and magnetization. We note that 
& = = in the limit of parallel propagating waves. It is therefore necessary to use 
other wave geometries in order to obtain an effect from the QED corrections. Note that this 
null-result should be expected on physical grounds, since successive Lorentz boosts along the 
direction of propagation would decrease the amplitude of all waves without limit. However, 
as shown in Refs. [3,4], a resonant nonlinear interaction between three waves is possible in 
bounded regions, as will be considered below. 

The most common approach [se e.g. Refs. [8]- [13]] to calculate the interaction strength 
is to use the general Lagrangian, apply the variational principle 5 f J&End 3 rdt = to find 
the influence of the QED terms in Maxwells equations (where the Lagrangian should be 
expressed in terms of the 4-potential), and proceed from there using standard techniques for 
weakly nonlinear waves. However, as shown in Ref . [4] , the algebra is significantly reduced if 
one make an ansatz for the potentials corresponding to the field geometry of consideration, 
and then derive the evolution equations directly from the variational principle, without using 
Maxwells equations as an intermediate step. 

Following the later approach we start by considering three interacting waves in a cavity 
with the shape of a rectangular prism. We let one of the corners lie in the origin, and we 
let the opposite corner have coordinates (xq, yo, z ). In practice we are interested in a shape 
where zq 3> xq, yo but this assumption will not be used in the calculations. We let the large 
amplitude pump modes have vector potentials of the form 

a a ■ ( nx \ ■ ( niTTz\ . . ^ . . 

Ai = Ai sin I — I sin I I exp(— \uj\t)y + c.c (2) 



x o J V z o 



. ny \ ( ri2ftz , 
A 2 = A 2 sin — sin I exp(— icj 2 £)x + c.c. (3) 



and 

A 2 = A 2 sin ( 

VoJ V zo 

where c.c. denotes complex conjugate, and we chose the radiation gauge $ = 0. It is easily 
checked that the corresponding fields 



/ 7T \ / itx\ ( n 1 / Kz\ 

B lz = — A 1 cos — sin exp(— iu;it) + c.c. , (4a) 

\x J \x J V z o J 

f n 1 -n-\ . { nx\ (niitz\ . . . 

n lx = — Ai sin — cos exp(— \uj-\t) + c.c. , (4b) 

V zo J \xo J V z J 

t-, a ( nx\ . ( uykz\ , . . . . 

Ely = lUiAi sin I — I sin I J exp(— lUit) + c.c. , (4c) 

together with = n\-K 2 c 2 / Zq + ir 2 c 2 /xl, and 
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_, / 7T \ I 7r y\ . / U 2 7TZ \ . . . . . 

i* 2z = - — U 2 cos — sin exp(-icu 2 t) + c.c. , (5a) 

\yoJ \yoJ \ z j 

( n 2 7r\ / 7iy\ ( n 2 nz\ 

B 2y = A 2 sm — cos exp(-iu; 2 t) + c.c. , (5b) 

V z o J \yoJ \ z o J 

7-. f 7[ y\ f n 2 irz\ 

E lx = \uj 2 A 2 sin — sin exp(— iu 2 t) + c.c. , (5c) 

\yoJ V -o J 

together with u 2 = n 2 n 2 c 2 / z 2 + n 2 c 2 jy\ are proper eigenmodes fulfilling Maxwells equations 
and the standard boundary conditions. Similarly we assume that the mode to be excited 
can be described by a vector potential 



A 3 = A 3 sin ( — ) sin ( — - — ) exp(— iu; 3 t)x + c.c. (6) 
V 2/o / V z o J 

where uo 2 = n 2 7i 2 c 2 /z 2 + 7c 2 c 2 /y 2 , in which case we get fields of the same form as in Eqs. 
(5a)- (5c). Since the QED terms are fourth order in the amplitude, the corresponding non- 
linearities are cubic, implying that u 3 = ±(2u 2 ± uji) or u 3 = ±(u 2 ± 2u>i) should hold for 
resonant interaction. We assume that the cavity dimensions are chosen such that a single 
eigenmode can be resonantly excited, and pick the alternative 

oo 3 = 2uoi - uo 2 (7) 

for definiteness. We note that when performing the variations 5A 3 , the lowest order terms 
proportional to vanish due to the dispersion relation, and we need to include terms 

due to the time dependence of the amplitude of the type A 3 d(8A%) / dt. For the fourth order 
QED corrections proportional to 5A* 3 , only terms proportional to A\A* 2 5A 3 survives the 
time integration, due to the frequency matching (7). After some algebra the corresponding 
evolution equation for mode 3 becomes: 

where the dimensionless coupling coefficient K rec is 
1 



K rec = — < -(-, +)2u 1 U 2 U 3 + 2nC 

n 2 n 3 



2^2 



.1 n~i 



0^0 



1 2 2 
+ 7T C 



3 n\n 2 n 3 n\ 

x oyo z o yo z o 
hn\uj2UJz 5u 2 n 2 n 3 



' 2u 2 uj 3 2u 2 \ 7ujiu 2 n 2 n 3 lu x u 3 n x n 2 

V ' ') 1 ,,2 ' „,2 I ' \'i ) ~2 ~r V ' J 2 



x o yo J Z Z 



(9) 



and we have added a phenomeno logical linear damping term represented by 7. When deriv- 
ing Eq. (9) we have assumed one of the following mode number matchings 

2ni - n 2 + n 3 = (10) 
2 ni +n 2 -n 3 = (11) 
2m - n 2 - n 3 = (12) 
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in order for the QED corrections terms to survive the z-integration. The three different sign 
alternatives in (9) correspond to the mode number matching options (10)-(12) respectively. 
Given experimental data for possible values of the pump field strengths and the damping 
coefficient inside cavities, the saturation level of the excited mode can be determined from 
Eq. (8). The possibilities for detection of photon-photon scattering using currently available 
performance on microwave cavities will be discussed in the final section of the paper. 



III. GRAVITATIONAL INTERACTION 

Preliminaries. In vacuum, a linearized gravitational wave can be transformed into the 
transverse and traceless (TT) gauge. Then we have the following line-element 

ds 2 = -c 2 dt 2 + dx 2 + [1 + dy 2 

+ [l-h + (0]dz 2 + 2h x (Odydz , (13) 

where £ = x — ct, c is the speed of light in vacuum, and \h+\, \h x \ <C 1. 

Neglecting terms proportional to derivatives of h + and h x (since the gravitational fre- 
quency is assumed small), the wave equation for the magnetic field is [5] 



1 d " ^2 
V 

c 2 dt 2 



B = 



. d 2 d 2 \ , d 2 

h + — - — + hy 



dy 2 dz 2 J dydz 



B, (14) 



and similarly for the electric field. 

In an isotropic dielectric medium with permittivity e different from the vacuum perme- 
ability, the equation (14) still holds, simply if we replace c by c/n = cj ^fe~ r in the above 
expressions, where £ — £ r £Q. For the moment, we will neglect mechanical effects, i.e., effects 
which are associated with the varying coordinate distance of the cavity walls due to the 
gravitational wave. 

Cavity design. The coupling of two electromagnetic modes and a gravitational wave in a 
cavity will depend strongly on the geometry of the electromagnetic eigenfunctions. It turns 
out that we can greatly magnify the coupling, as compared to a rectangular prism geometry, 
by varying the cross-section of the cavity, or by filling the cavity partially with a dielectric 
medium. The former case is of more interest from a practical point of view, since a vacuum 
cavity implies better detector performance, but we will here consider the later case since it 
can be handled analytically. 

Specifically, we choose a rectangular cross-section (side lengths Xq and yo) , and we divide 
the length of the cavity into three regions. Region 1 has length li (occupying the region 
—l\ < z < 0) and a refractive index ri\. Region 2 has length l 2 (occupying the region 
< z < l 2 ), with a refractive index n 2 , while region 3 consists of vacuum and has length Z 3 
(occupying the region l 2 < z < l 3 + l 2 ). We will also use I — h + l 2 + Z 3 for the total length. 
The cavity is supposed to have positive coordinates, with one of the corners coinciding with 
the origin. Furthermore, we require that l 2 li, and that the wave number in region 2 is 
less than in region 1. The reason for this arrangement is twofold. Firstly, we want to obtain 
a large coupling between the wave modes, and secondly we want an efficient filtering of the 
eigenmode with the lower frequency in region three. 
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The first step is to analyze the linear eigenmodes in this system. Those with the lowest 
frequencies are modes of the type 



_ icuxq ~ / mux \ . _ • . 

E y = B Z j sin I J sm[kjZ + (pj\e , 

17171 V Xq I 

5 2 = 5^- cos ' » ~~r;„ ~ , ,„ 



Z? 7 



T17TX \ 

B? A sin ( ) cos[/cj2; + v^Je" 



■^B zj sin — r '- ~ 1 

m7r \ Xq / 



(15a) 
(15b) 
(15c) 



in regions j — 1, 2 and 3, respectively, where the wave in region 3 is a standing wave. 
Furthermore, in region 3 we may also have a decaying wave 



E„ 



±? z3 sin ( I sinh[/c 3 2 + y? 3 Je , 



77T.7T 

B z = B z3 cos 



m7rx\ 

smh[k 



—iujt 



Xq J 

D 5 . ( mirx\ 

B x = B z3 sin 

m7r \ x / 



3^ + <Ps]e 



coshfA^z + ips}e' 



-iuit 



(16a) 
(16b) 
(16c) 



Using standard boundary conditions, it is straightforward to perform most of the eigenmode 
calculations analytically. Once the wavenumbers are calculated for an eigenmode, the re- 
lation between the amplitudes in the three regions is found, and thereby the mode profile. 
We are specifically interested in the shift from decaying to oscillatory behavior in region 3, 
and we denote the highest frequency which is decaying in region 3 with index a, and the 
wave number and decay coefficient with k la , k 2a and k 3a respectively. Similarly, the next 
frequency, which is oscillatory in both regions, is denoted by index b. If we have / ^> x 
(and m the same) these two frequencies will be very close, and a gravitational wave which 
has a frequency equal to the difference between the electromagnetic modes causes a small 
coupling between these modes. An example of two such eigenmodes is shown in fig. 1. 




350 



FIG. 1. Mode profiles for the eigenmodes are shown for the parameter values rt\ 
n 2 = 1.0005, hm/x0 = 143, l 2 m/x0 = 12.3 and l 3 m/x0 = 184. 



1.22, 



-Bnorma is the solid line and 



Bnormb is the dotted line, and I c = 0.27. 
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We define the eigenmodes to have the form B zab = B ab (t)B noTmab (z, y, z), where B a b is a 
time-dependent amplitude and the normalized eigenmodes -B n orma,fe fulfill f v |-B n orma,b| 2 dV = 
V. We let all electromagnetic field components be of the form A = A a (r) exp(— iuo a t) + 
Ab(r) exp(-iubt) + c.c, where c.c. stands for complex conjugate, and the indices stand for 
the eigenmodes discussed above. The gravitational perturbation can be approximated by 
^+,x = h +jX exp(-iujgt) + c.c, where we neglect the spatial dependency of h +;X altogether, 
since the gravitational wavelength is assumed to be much longer than all of the cavity 
dimensions. If we consider a binary system of two black holes fairly close to collapse, the 
gravitational frequency will not be an exact constant, but will increase slowly. During a 
certain interval in time, the frequency matching condition 

uj b = u a + u g , (17) 

will be approximately fulfilled. Given the wave equation (14), and the above ansatz we find 
after integrating over the length of the cavity 

' 2 '^ ( d ~ l) B b = -h + kl a B a I c , (18) 



c 2 \dt 
where 

'-Fk/^'" (19) 

and we have added a phenomenological linear damping term represented by 7. Thus we note 
that for the given geometry, only the /i + -polarization gives a mode-coupling. (Rotating the 
cavity 7r/4 around the x-axis will instead give coupling to the /i x -polarization.) Furthermore, 
if we consider propagation in a different angle to the cavity , the result will be slightly 
modified. Calculations of the eigen-mode parameters show that I c may be different from 
zero when ni 2 7^ 1, and generally I c of the order of unity can be obtained, see fig. 1 for 
an example. From Eq. (18), we find that the saturated value of the gravitationally excited 
mode is _ 

~ h + k la c B a 

-Dfeat - ~ i-c ■ (/UJ 

27^6 

In fig. 1 it is shown that we can get an appreciable mode-coupling constant I c for a cavity 
filled with materials with different dielectric constants, and it is of much interest whether the 
same can be achieved in a pure vacuum cavity. As seen by Eq. (19), the coupling is essentially 
determined by the wave numbers of the modes, given by k = (m 2 n 2 /xQ — n 2 u; 2 /c 2 ) 1//2 . Thus 
by adjusting the width xq in a vacuum cavity, we may get the same variations in the wave 
numbers as when varying the index of refraction n. The translation of our results to a 
vacuum cavity with a varying width is not completely accurate, however. Firstly, when 
varying x , the mode-dependence on x and z does not exactly factorize, in particular close 
to the change in width. Secondly, the contribution to the coupling I c in each section becomes 
proportional to the corresponding volume, and thereby also to the cross-section. However, 
since most of the contribution to the integral in Eq. (19) comes from region 1, our results 
can still be approximately translated to the case of a vacuum cavity, by varying x instead 
of n such as to get the same wavenumber as in our above example. Thus we conclude that 
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our discussion of the sensitivity of a cavity based detector given below, can be based on a 
vacuum cavity rather than a "dielectric cavity", where the former case is preferred due to 
the much smaller dissipation level. 



IV. DISCUSSION OF THE DETECTION SENSITIVITY 

For both cases considered above, as can be seen from equations (8) and (20) respectively, 
the quality factor Q of the cavity (i.e. the damping time/periodtime) and the allowed field 
strength are the main parameters that determine the saturation level of the excited mode. 
For superconducting niobium cavities large surface field strengths (E ~ 60 MV/m) and high 
quality factors Q ~ 10 10 can be obtained simultaneously [14]. For the case of photon-photon 
scattering this typically corresponds to a few excited microwave photons in the new mode. 
Provided the number of excited photons beats the thermal noise level (N th ~ kT/hu 3 ) a few 
microwave photons is enough for detection, see e.g. [15] . However, detection of such a small 
signal can only take place if the pump waves are filtered out. This can be done directly in 
the cavity using a "filtering geometry", i.e. a cavity with a variable cross-section, where only 
the excited mode has a frequency above cut-off in the detection region of the cavity. Our 
conclusion is that detection of photon-photon scattering is possible in a microwave cavity 
using current technology, provided the equipment has state of the art performance, and a 
filtering geometry is applied. 

The considerations for the gravitational wave detector are similar to the photon-photon 
scattering case. However, we must also take the length variations of the detector due to 
acoustic thermal noise into account. As a consequence, an efficient detector must then have 
a rather high mechanical quality factor Q mcc ~ 6 x 10 6 (where Q m ec describes the relative 
damping of acoustic oscillations), for sensitive detection to be possible. Furthermore, we 
must also be aware that the gravitational wave sources of most interest - binary systems 
close to collapse - produces radiation with a finite frequency chirp, implying a finite time 
of coherent interaction, typically fractions of a second. If the above aspects are taken into 
account, we find that a cavity based detector of a few meters length can detect metric 
perturbations of the order h min ~ 7 x 10~ 23 (See Ref. [5] for a more thorough discussion of 
the detection level.). 
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